Spinwaves and large electron-phonon coupling near the metal-insulator transition in 

hole-doped high-Tc oxides. 



DPMC, University of Geneva, 24 



o 
o 

(N 

X> 

<D 

O 

(N 

O 

o 

u • 

Of 

3 ; 

czj ; 



^3 

G 

O 

o 



> 
m 

(N 
O 



c3 



i 

c 

o 
o 



X 



T. Jarlborg 

ljuai Ernest- Ansermet, CH-1211 Geneve 4, Switzerland 
(February 1, 2008) 



A theory for spin and lattice couplings to the electronic states in high-Tc oxides is presented, with 
HgBa2Cu04 as an example. A simple analytical model suggests that the barrel band is sensitive to 
potential perturbations with long-range Fourier components, and that gaps can be opened at the 
Fermi energy. This is confirmed in self-consistent band calculations for modulations in elongated 
supercells, where spin waves and phonon coupling compete for equal q- vectors. It is argued that 
the wavelengths of spin-wave modulations and of the phonon modes with large coupling, depend 
on the doping. This mechanism supports the idea that pseudogaps are caused by stripe like spin- 
modulations in underdoped systems, while superconductivity is attributed to enhanced electron- 
phonon coupling for long wavelength phonons. 



The mechanism behind high-Tc superconductivity in 
copper oxides is still unknown. From the vast amount of 
experimental information that has been collected from 
high-Tc oxides, one can recall some facts that seem 
to distinguish Uupfe materials from ordinary low-Tc 
superconductorsEra Phase sensitive tests of the pairing 
indicate d-wave symmetry of the superconducting gapQ. 
This is also implied from photoemission and tunneling 
measurements where tha-gap is highly anisotropic over 
the Fermi Surface (FS)ErEI. Such measurements detect 
pseudogaps, which extend well above Tc in underdoped 
systems. They appear to have the same symmetry as the 
superconducting gap, making the FS visible only along 
the diagonal. The isotope effect, which is small and vari- 
able, exists also for the pseudogapfl. Spin fluctuations 
follow the same type of antiferromagnetic (AFM) orienta- 
tion of ordered Cu moments that can be seen in undoped, 
insulating copper oxides. Separation of metallic and inst 
lating regions form static stripes for some compositions 
but they might be dynamic in other cases. Band the- 
ory based on the density functional approximation fails 
to describe the undoped, insulating cuprates, while the 
bands and FS of metallic, superconducting cuprates agree 
well with experimental data. The elcclmn- phonon cou- 
pling (A) calculated from band theorytlE], is quite large 
in view of the low density-of-states (DOS) in these mate- 
rials, but still not high enough to explain the high Tc's. 
Here we perform self-consistent electronic structure cal- 
culations based on the local (spin) density approximation 
(LDA)Eil, for long wavelength excitations of either spin 
fluctuations or phonons, by use of large supercells. As 
will be shown, spin waves and phonons seem to com- 
pete for the same q-vectors, and the perturbations will 
be large at the FS, leading to a tendency for gap forma- 
tion in both cases. We focus on HgBa 2 Cu04 (HBCO), 
since it is relatively simple with one CuO-plane and only 
one cylindrical FS, a so-called 'barrel'. Such a piece of FS 
is common to all high-Tc oxides and probably involved 
in the mechanism for superconductivity. The lattice con- 
stant a is 7.32 a.u.. The bands are calculated using the 
linear muffin-tin orbital method, where 5 empty spheres 



are introduced to account for non-spherical potential cor- 
rections in the most open parts of the structure. The 
structure, details of the method ancUaaHd results for the 
basic unit cell, are described earlierli3il3. Only Cu, and 
to some extent the planar O, have a significant DOS at 
Ep. The band structure of the normal cell is shown in 
fig. |l| for k-points in the k z =0 plane. The FS is a cylin- 
drical 'barrel' centered at the M-point. The FS reaches 
the zone boundary at X for a doping of about 0.25 hole 
per formula unit (f.u.). The cell of the AFM insulating 
case, has two Cu in the xy-plane with opposite magnetic 
moments. This cell of two f.u. has 26 sites (16 atoms and 
10 empty spheres). The AFM state can be reached by 
introducing large correlation, increased exchange or by 
application of external field. Here we include a staggered 
field (positive and negative with maximum amplitude of 
30 mRy) on the two Cu in the AFM cell. The undoped 
material is insulating with this field, with a gap of 10 
mRy and local moments of about 0.5 \ib- The critical 
field for having a gap is about 20 mRy. A uniform field 
within each Cu-sphere makes the entire Cu-bands to be 
shifted, so that part of the upper DOS-edge of the Cu-d 
band will get closer to Ep with increasing field. 

The 26-site cells are put together in the [1, In- 
direction, to form larger cells in order to study q-waves 
modulated along the diagonal (with respect to the Cu- 
O bonding). Four supercells containing 104, 130, 156 
and 208 sites are considered, corresponding to a maxi- 
mal wavelength of 8y2a. For modulations along [1,0,0] 
it necessary to have a two-layer wide unit cell (along y) 
because of the underlying AFM order. A cell of 208 sites 
will describe spin modulation with wavelength 8a along 
x. For studies of phonon modulations without consider- 
ation of spins, we can use narrow cells (width a along 
y), so that a 104-site cell can describe a modulation with 
period [8a, 0,0]. In each case we calculate the bands in 
24 k-points in the respective irreducible Brillouin Zones 
(BZ). The DOS is determined by from tetrahedron k- 
point integration, where it is assumed that the energies 
vary linearly within each tetrahedron. The use of few 
k-points may introduce artifacts in the DOS, especially 
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near gaps where the bands deviate from linearity. 

The z-dispersion of the barrel band is very small, but 
the dispersion is large in the xy-plane, as along the diag- 
onal r — M (cf. fig. A modulation of the potential 
V q - sin(gr), expressed by the Fourier coefficient V q , where 
q is along [1,1,0], will open a gap for a free electron like 
band. From the well-known solution for bands near the 
zone boundary in semiconductorsE3, 

E(k) = \{e k + e k - q ) ± \^{e k - e fe _,) 2 + 4V? (1) 

which at the zone boundary k — ^q, becomes E(q) = 
ei q ± Vq. Thus, there are two states with a gap, 2V q , in 

between, at the new BZ-limit at \q. The square of the 
wave function belonging to the state below the gap, has 
a modulation in phase with the potential modulation so 
that its amplitude is largest where the potential is most 
attractive. The square of the wave function for the state 
above the gap has the same modulation, but its phase 
is shifted, so that the large amplitude coincides with the 
most repulsive part of the potential. 

Consider first AFM cases. The AFM state in the un- 
doped case is partly stabilized because of the opening of 
the gap at Ep. Hole doping made in a rigid-band man- 
ner on the same DOS lead to metallic states, without the 
stabilizing effect of the gap. As will be seen, modula- 
tions of the AFM moments will open gaps at different 
doping levels, with possible energy gains. An example 
of a modulation is when the Cu sites projected on the 
[1,1,0] direction have moments; zero, up, down, up, zero, 
down, up, down. Such a "stripe" has ±V q components 
in the spin-potentials, with a wavelength covering totally 
8 Cu layers. If Ep is far below (or far above) the gap, 
both states will be occupied (or unoccupied) with no gain 
in total energy. It is only if Ep is near the gap that 
there will be a possibility to gain total energy. The oc- 
cupied state reduce its kinetic energy, while the raise in 
energy of the upper, empty state does not contribute to 
the total energy. The process is self-supporting via the 
exchange potential, since the increased spin density in re- 
gions where the potential is spin splitted will increase the 
potential splitting further. Thus, the argument is that 
the CuO-plane systems can chose the wavelength of the 
modulation according to its doping to make an optimal 
gain in total energy. This mechanism always opens the 
gap (or pseudogap) close to Ep, leading to a correlation 
between the wavelength and doping. These arguments 
can be used as long as there are no other bands near Ep, 
because they would normally involve an increase of the 
total energy when a modulation is imposed. 

Similar arguments can be used to discuss the coupling 
between phonon modes and the electronic states near Ep. 
For example, a phonon induced modulation of the poten- 
tial on Cu sites will also have its Fourier component V q . 
(Since the local Cu-DOS is the highest near Ep, it will be 
largest effects from potential changes on Cu-sites). The 
(non spin-polarized) bands will open a gap as for each of 



the spin bands in the AFM case, if the gap is close to 
Ep. As before, only one particular wavelength will be 
concerned, depending on the doping. It is worth noting 
that for a given doping, the same q- vector is involved for 
both spin waves and phonons. This mechanism will be 
dynamic for phonons, unless the gain in energy will make 
the phonon completely soft to make a lattice instability. 
Spin waves could be pinned to impurities and therefore 
be static, but coupling between spin waves and phonons 
are likely since they compete for the same q- vectors. 

Next, we consider the band results. The basic results 
from the simple, one-dimensional model, are confirmed 
by the band calculations, and the gap can be formed over 
the whole FS for varying doping. For the cases with AFM 
modulations, we first apply sinus-like modulations in the 
diagonal direction. The staggered fields (on Cu sites) 
have the same maximal amplitude as for the undoped in- 
sulating case. The undoped system contains 63 electrons 
per f.u.. However, the gaps (or pseudogaps) found in the 
four cases with elongated cells are always found at one 
state (1 electron of each spin) below this electron filling. 
This means at band fillings of 502, 628, 754 and 1006 
electrons, for the four supercells of increasing lengths. 
The doping is adjusted to coincide with these fillings, by 
use of the virtual crystal approximation. The DOS near 
Ep with and without AFM modulation, are shown in fig. 
||for the 104-site cell. The gap is small, like a pseudogap, 
and it coincides with Ep when the doping is 0.25 holes 
per f.u., corresponding to the particular q. Modulations 
in the longer cells give similar small gaps, or pseudogaps, 
but at positions corresponding to decreased hole doping, 
at 0.2, 0.1667 and 0.125 holes per f.u., for the cells with 
130, 156 and 208 sites, respectively. The opening of the 
pseudogap starts along the diagonal T — M . Larger fields 
are needed to form a complete gap extending over the 
whole FS. Thus, in agreement with the simple model, 
the band results show that the barrel band can be cut to 
form gaps at different positions along Y — M, depending 
on the Fourier component of the potential perturbation. 

AFM modulations in the x-direction is made for the 
208-site cell, where the pseudogap is found for a dop- 
ing of 0.125 hole per f.u., as for the longest cell along 
the diagonal direction. By folding the k-points in the 
BZ back to the original BZ, it is seen that the partial 
opening of a gap is largest far from the diagonal. This 
is consistent with the observations made by photoemis- 
sion, where the FS seems to "disappear" except near the 
diagonalfl. (It,is assumed that stripes are running both 
along x and yfiJ, making their effect superimposed.) Also 
the observation of static stripes along [1,0,0] in one high- 
T c materials suggests that [l,0,0]-modulations are more 
probable than diagonal ones. This can be understood 
from the band results as follows: The larger z-dispersion 
near off-diagonal points on the FS can be translated into 
higher DOS on the FS for these points, so that a band 
opening at those points will lead to a larger gain in total 
energy than if the gap opened along the diagonal. As the 
orientation in real space can be associated with the dif- 
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ferent FS points, it is therefore preferable for the system 
to chose the [l,0,0]-orientation leading to the largest gain 
in energy. 

These spin-polarized band calculations are not ab- 
initio, and they cannot determine the precise form and 
amplitude of the modulations in the doped cases. One 
can expect that some kind of modulation is present, 
which in this band description produces gaps and stripes. 
Other theories based on strong on-site correlation have 
been showa to lead to stripe formation and spin-charge 
separationEl. Here, in a band description, the spin, charge 
and local DOS information show that the charges on Cu 
of different polarization are very similar. The DOS at 
Ep is largest on the non-polarized Cu-sites, but still siz- 
able on the polarized Cu sites, at least for one of the 
spins. Thus, conduction is not only restricted to the zero- 
moment part of such stripes. 

In the present model, there is a close connection be- 
tween the doping and the wavelength of the stripe. The 
lower the hole doping, the longer wavelengths. Further, 
it is probable that long wavelength modulations are more 
stable than short ones. This is because at least two Cu 
sites have no moment in each modulation, i.e. no ex- 
change energy can be gained from a polarization of these 
two sites. The fraction of zero moment Cu sites is larger 
for short wavelength modulations, making them less sta- 
ble. This cannot be verified by ab-initio total energy 
calculations, but for the short modulation of only six 
Cu-layers along [1,0,0] it is more difficult to open a gap 
even for large field amplitudes. Hybridization between 
near neighbors in combination with rapid variations of 
the potential over only a few layers will prevent an easy 
opening of a gap. In addition there is a topologic change 
of the FS at the X-point when the doping approaches 
0.25. Therefore, a possibility is that pseudogaps, caused 
by spin waves, cannot exist for too large doping. This is 
consistent with observations of pseudogaps only in under- 
doped high-Tc oxides. A. critical concentration of about 
0.19 has been suggested!! This is to be compared with 
the values 0.125 for the shortest modulation along [1,0,0] 
and 0.25 along [1,1,0] to produce a "good" pseudogap in 
the calculations. 

The non-spinpolarized band calculations for phonon 
modes are ab-initio within LDA, for the imposed lattice 
distortions. Coupling between spin waves and phonons 
is plausible, but we do not study such cases here. Dif- 
ferent types of lattice distortions are considered, mod- 
ulated along [1,0,0] and [1,1,0] in the elongated super- 
cells. Several distortions have a direct influence on the 
Cu potentials like the change the distance between Cu 
and apical-O, or between Ba and the CuO plane. The 
strongest effects occur for "breathing" modes, where the 
four planar O sites are moved alternatively towards or 
away from the central Cu. The modulations are cho- 
sen to be sinus-like. Distortion amplitudes of the order 
0.01-a are expected at low temperature for typical force 
constants and zero-point motionsL3, but in order to en- 
hance the effects to be seen in the DOS, we also consider 



amplitudes of 0.03-a. 

The DOS for the 104 site cell is shown in fig. ||. For 
smaller distortion amplitudes the effect is smaller and 
hardly seen on the DOS. By comparing the back-folded 
k-points with the points in the normal BZ, one can verify 
that the effects due to modulations along the diagonal 
or along the CuO-bonding, indeed are concentrated to 
the expected parts of the FS. Moreover, if the sinus-like 
modulation is replaced by the absolute value of sinus, one 
finds about one order of magnitude smaller energy shifts 
near Ep. Results for the other cells are similar, so that 
the distorted structures show some different DOS com- 
pared to the undistorted cases, but no real gaps are de- 
veloped. However, strong changes of the energies means 
that A is large for these modes. Full calculations of A 
will not be presented here, but a few observations can 
be made. As was argued above, the wavelengths are de- 
termined via doping so that only energies close to Ep 
are changing. Calculations of electEoa-phonon coupling 
made for the undoped compositionsBE3 will miss this ef- 
fect, since the sensitivity to V q of a certain wavelength 
appears below the position of Ep for the undoped sys- 
tem. (The wavelength corresponding to no doping would 
be infinite.) The largest change in energy from the model 
above, is equal to V q , i.e. equal to the maximal ampli- 
tude of a sinusoidal potential perturbation. Two exam- 
ples of O- and Ba-distortions with V q from the potential 
in the band calculation, give A in.-the range 1-2. This 
is from an approximate expressionlilJ for the long-range 
coupling, A = NV^/Y,K{^R) 2 , where N is the DOS 
at Ep per cell, if is a force constant0, AR the dis- 
placement of each site, so that K times the sum over 
all sites gives the total elastic energy of the deformation. 
Electronic screening in a good metal tend to diminish the 
long-range part of A, since charge relaxation close to a lo- 
cal lattice distortion will prevent any far-reaching change 
in Coulomb potentialt-i As these materials are close to 
a metal-insulator transition, it is possible that the mech- 
anism of screening is reduced. The mechanism should 
retain the property of forward scattering of A which is 
favorable to d-wave pairingEj. Strong coupling for large 
wavelengths will renormalize phonons near the zone cen- 
ter. It is not yet clear if low or high frequency phonons 
will contribute most to superconductivity. 

In conclusion, the simple model of gap formation, sup- 
ported by band results for HgBa2Cu04, suggests a con- 
nection between spin waves, stripes, pseudogaps and 
strong electron-phonon coupling. The reason is that the 
barrel band is sensitive to short-q perturbations, V q . This 
scenario explains a number of features of the high-Tc ox- 
ides mentioned above. The origin of the pseudogaps is the 
stripe like spin-waves, which are stable only below a cer- 
tain doping, which in our model means above a certain 
wavelength of the spin modulation. Whether they are 
dynamic spin-fluctuations, static or coupled to phonons, 
cannot yet be determined from the present calculations, 
which are static. Energy considerations suggest that the 
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spin waves are more probable at long wavelengths, i.e. in 
undcrdoped cases. Decreasing stability of the spin waves 
at short wavelength, i.e. for increasing hole doping, sug- 
gests that they and their role diminish with increasing 
doping. If this picture is true, and spinfluctuations are 
harmful to superconductivity, one can expect that by 
applying pressure to samples with little doping it will 
be possible to suppress the fluctuations and reach higher 
Tc- The effect of pressure should be smaller in overdoped 
samples. The band is equally sensitive to perturbations 
V q , produced by phonons. The amplitudes of V q due to 
phonons of realistic deformation amplitudes at low T, are 
smaller than in the case of AFM modulation, and fully 
developed gaps are normally not visible in the DOS. But 
a particular V q induced by phonons favors band shifts at 
some sections of the FS and a large A. Any perturbation, 
phonons, spin or valence fluctuations leading to a similar 
V q component will show similar effects on the electronic 
states. Competition between different mechanisms can 
make interpretation of experiments difficult. 
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FIG. 1. Band structure of HgBa2Cu04 along symmetry 
directions in the k z =0 plane of the tetragonal BZ. The Fermi 
energy is indicated by the broken line. 



J. Millis, Science 288, 468, (2000), and 
Prog. Phys. 62, 61, (1999), 
Physica C 349, 53 (2001), 
Rev. Lett. 85, 182, 



1 J. Orenstein and A. 
references therein. 

2 T.Timusk and B.Statt, Rep. 
and references therein. 

3 J.L. Tallon and J.W. Loram, 
and references therein. 

4 C.C. Tsuei and J.R. Kirtley, Phys. 
(2000). 

5 M.R. Norman et al, Nature 392, 157, (1998). 

6 Ch. Renner et al, Phys. Rev. Lett 80, 149, (1998). 

7 J. Hofer et al, Phys. Rev. Lett 84, 4192, (2000). 

8 J. Tranquada et al, Phys. Rev. Lett 78, 338, (1997). 

9 H. Krakauer et al, Phys. Rev. B 47, 1002, (1993). 

10 T. Jarlborg and G. Santi, Physica C 329, 243 (2000), and 
references therein. 

11 W. Kohn and L.J. Sham, Phys. Rev. 140, A1133, (1965); 
O. Gunnarsson and B.I Lundquist, Phys. Rev. B13, 4274, 
(1976). 

12 B. Barbiellini and T. Jarlborg, Phys. Rev. B50, 3239 
(1994). 

13 J.M. Ziman, "Principles of the Theory of Solids", Cam- 
bridge (1971), p71. 

14 T. Jarlborg, Phys. Rev. B 61, 9887, (2000). 

15 G. Santi et al, J. Supercond. 8, 405 (1995), G. Santi and 
T. Jarlborg, Int. J. Mod. Phys. 13, 3566, (1999). 



104-site cell along [1,1,0] 



normal (para) 

with phonon distortion (para) 

modulated AFM (spin-pol) 




(E-E p ) (Ry) 



FIG. 2. Paramagnetic DOS for a 104-site unitcell of 
HBCO (doped with 0.25 holes per f.u.) along [1,1,0] (thin 
line), paramagnetic DOS when the structure contains a sinus 
modulation of Ba z positions (heavy line) and the DOS for a 
sinus-modulated AFM state along the cell (broken line) . The 
latter is obtained via applied fields (maximum amplitude 30 
mRy) on the Cu sites. A higher field will open the gap com- 
pletely, but also move the Cu-d band edge closer to Ef- 
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FIG. 3. Paramagnetic (thin line) and AFM (heavy line) 
DOS for a 208-site unitcell of HBCO (doped with 0.125 holes 
per f.u.) along [1,0,0]. 



